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Abstract 
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Submitted to JHEP 



^Research supported by the US Department of Energy under grant DE-FG03-91ER40662. 
''^Research supported by the US Department of Energy under grant DE-AC03-76SF00515. 
''Laboratory of the Direction des Sciences de la Matiere of the Commissariat a, I'Energie Atomique of France. 



1 Introduction 



The quest for more precise theoretical predictions for particle production at colliders requires cal- 
culations to ever-higher orders in the perturbative expansion of field theories. At tree level and at 
one loop, workers have made substantial progress in recent years in computing multi-parton QCD 
scattering amplitudes ||. Several important quantities, such as the total cross section for e~^e~ 
annihilation into hadrons and the QCD /3-function, have been calculated up to four loops j^. In 
contrast, essentially no higher- loop fixed-order results are available for multi-parton processes de- 
pending on more than one kinematic variable. In particular, the two-loop amplitudes necessary 
for reducing the theoretical uncertainty in the measurement of Og from e~^e~ jet rates and other 
event-shape variables remain uncalculated. 

There has, however, been some progress in developing general formalisms for the computation 
of next-to-next-to-leading-order (NNLO) jet cross sections in QCD. This includes the structure of 
the infrared singularities that arise when partons are unresolved and the explicit evaluation of 
two-loop integrals [||, ^ associated with massless processes. In the special case of maximal = 4 
supersymmetry, one of the authors and his collaborators have computed the two-loop four-point 
amplitudes in terms of scalar double box integrals [0]. 

In this paper we compute a non-supersymmetric two-loop QCD amplitude. Our method for 
obtaining the amplitude is based on using unitarity to determine its functional form. This program 
has been carried out for many one-loop and a small of number of two-loop amplitudes 0, ^ p^] . 
The latter computations were for the special case of theories with A = 4 or A = 8 supersymmetry. 
The two-loop amplitude considered here is more difficult to obtain because one cannot directly rely 
on supersymmetry. We will have to evaluate the unitarity cuts with more care, by working in 
dimensional regularization with arbitrary dimension D = 4 — 2e, taking the limit e ^ only at the 
end. 

The particular process we study is g~g~ g~^g~^ in pure gauge theory, where g~ {g~^) denotes a 
negative (positive) helicity gluon. In a helicity convention in which all gluons are treated as outgoing, 
the amplitude is that for four identical-helicity gluons, — > g~^g^g^g^- This helicity amplitude 
vanishes at tree level [[Til], and consequently the two- loop term we compute here contributes to the 
cross section only at one order beyond NNLO. Nevertheless, the techniques used here should be 
applicable to the helicity configurations that do contribute at NNLO. 

This paper is organized as follows. In section |2| we give an overview of the unitarity technique 
as applied to two- loop amplitudes. Section ^ describes the construction of some useful auxiliary 
amplitudes where scalars replace some of the gluons circulating in the loops. The full two-loop 
four-gluon amplitude is presented in section ^. In section ^ we verify that this amplitude has the 
combined ultraviolet and infrared divergence structure expected on general grounds [12, An 
appendix contains the expansion in e, through 0{e^), of the two-loop integrals appearing in the 
amplitude, expressed in terms of standard functions (poly logarithms) . 
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2 Overview of method 



We reconstruct the two-loop four-gluon amplitude from its unitarity cuts, by evaluating the latter 
in dimensional regularization |14] with D = 4 — 2e. Traditional applications of unitarity, via dis- 



persion relations, suffer in general from subtraction ambiguities. These ambiguities are related to 
the appearance of rational functions with vanishing imaginary parts, R{Si), where Si = {s, t,u, . . .} 
are the kinematic variables for the amplitude. However, dimensionally-regulated amplitudes for 
massless particles necessarily acquire a factor of (— Sj)"^ for each loop, from the measure J d^L and 
dimensional analysis. For small e, {—S.i)~'^ R{Si) = R{Si) — e ln{— Si) R{Si) + ©(e^), so every term 
has an imaginary part (for some Si > 0), though not necessarily in those terms which survive as 
e ^ 0. Thus, the unitarity cuts evaluated to 0(e) provide sufficient information for the complete 
reconstruction of an amplitude through 0{e^), subject only to the usual prescription dependence as- 
sociated with renormalization |^, ^. The subtraction ambiguities that arise in traditional dispersion 
relations are related to the non-convergence of dispersion integrals. A dimensional regulator makes 
such integrals well-defined and correspondingly eliminates the subtraction ambiguities. In a sense, 
we use dimensional regularization as a calculational tool, beyond its usual role as an infrared and 
ultraviolet regulator. 

We make use of the fact that the amplitudes we are computing are also expressible in terms 
of Feynman diagrams. They can therefore be expressed as a sum of loop integrals multiplied by 
kinematic coefficients. We can therefore bypass the dispersion integrals by identifying the loop 
integrals whose cuts match the appropriate products of tree or one-loop amplitudes. Our main 
task will then be to compute the integrands of the two-loop integrals from the tree or one-loop 
amplitudes. We obtain the full amplitude in terms of various loop integrals, by requiring consistency 
of the different cuts. The consistency conditions are relatively simple to implement because they 
hold before integrating over intermediate-state momenta. Because no integration is required, we 
can verify the resulting form numerically (to high accuracy), even when analytic simplification of 
some of the cuts is impractical. This approach is powerful because amplitudes on both sides of the 
cuts can be simplified, including cancellations dictated by gauge invariance, before their products 
are computed. This reduces greatly the number of terms in intermediate steps. The final results 
are quite compact when expressed in terms of loop-momentum integrals, thus simplifying the task 
of performing the integrals explicitly. The cut-based method thereby provides a short-cut to final 
expressions for amplitudes. 

Let us begin by considering the two-particle cuts of the full amplitude, including color factors. 
In the s channel, in general, there are two contributions to such cuts, as displayed in fig. [l|. These 
are, 



^f'°°P(l,2,3,4) 

J2 [AT%l,2,-i2,-h) X ^l-'°°P(^i,£2,3,4) +4-'°°P(l,2,-^2,-^i) X ^f^(4,^2,3,4) 



2-cut 



physical 
states 



(2.1) 

where the notation '2-cut' means cutting the lines corresponding to £±^2, taking the absorptive part. 
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and extracting the integrand of the resulting phase-space integral. On the right-hand side the legs 
labeled £1^2 are thus on shell in D dimensions, £f = 0. Here A^^^^"^, and Af^^°°^ are full n-point 

amplitudes, including all color factors. For each external leg we have abbreviated the dependence 
on the outgoing external momenta ki, color indices a^, and polarizations £i, by the label i. The 
sum over physical states crossing the cut contains a sum over colors. If the cut leg is a gluon, its 
associated polarization vector ee- must be transverse, ££. ■ ii = 0. In the sum over states, this can be 
imposed by a transverse projector Pf^i, acting on the Lorentz indices crossing the cut. 

We must also consider three-particle cuts. In the s-channel this cut is displayed in fig. and is 
given in terms of five-point tree amplitudes by, 



^^-^°°P(1,2,3,4) 



3- cut 



E -4r(i,2, 



(2.2) 



physical 

states 



where the notation '3-cut' means an operation analogous to that denoted by '2-cut', but cutting the 
three lines corresponding to ^1.2,3- On the right-hand side the corresponding legs are on shell. The 
two-loop amplitude will satisfy both eqs. ( ^.1| ) and ( |2.2[) . These s-channel cuts, together with the 
similar t- and it-channel ones, provide the analytic information necessary to obtain the complete four- 



point two- loop amplitudes. (The invariants s,t,u are the usual Mandelstam variables, s 
t = {ki + k^)"^, and u = {ki + ks)"^.) 



{ki + k2 





Figure 1: The s-channel two-particle cuts of a two-loop amplitude as products of tree and one-loop ampli- 
tudes. We take all external momenta to be outgoing. The dashed lines represent the cuts. 




Figure 2: The three-particle cut of a two-loop amplitude. 

The unitarity-based technique may also be viewed as an alternate way of evaluating sets of 
ordinary Feynman diagrams. It does this by collecting together gauge- invariant sets of terms which 
correspond to the residues of poles in the integrands. The poles are those of the propagators of the 
cut lines. This corresponds to a region of loop-momentum integration where the cut loop momenta 
go on shell and the corresponding internal lines become the intermediate states in a unitarity relation. 
From this point of view, we may consider even more restricted regions of loop momentum integration, 
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where additional internal lines go on shell (and, if they are gluons, become transverse as well). This 
amounts to imposing cut conditions on additional internal lines. 

In this vein it is useful for us to define a 'double' two-particle generalized cut for a two-loop four- 
point amplitude. This quantity, illustrated in fig. |^, is written in terms of on-shell tree amplitudes 
as 



1,2,3,4) 



2x2-cut 



'1,^2, 



physical 
states 



(2.3) 

where the on-shell conditions are again imposed on the ii, i = 1, 2, 3, 4 appearing on the right-hand 
side. This equation should not be interpreted as trying to take 'the imaginary part of an imaginary 
part'. Rather it should be understood in the sense of the previous paragraph as supplying information 
about the integrand of the two-loop amplitude. It supplies only part of the information contained 
in the usual two-particle cut ( p.ip , which effectively imposes only two kinematic constraints on the 
intermediate lines. However, it is simpler to evaluate because it is composed only of tree amplitudes. 
Its computation is thus a natural first step. As we shall see this double cut contains 'most' of the 
final result for the amplitude we are computing in this paper. 







E 









Figure 3: The s-channel double two-particle cut of a two-loop amplitude separates it into a product of three 
tree amplitudes. The dashed lines represent the generalized cuts. 



Consider the two- loop four-point amplitude with all particles in the adjoint representation of the 
gauge group SU{Nc). We may write its color decomposition as follows 

^4"'°°P(1,2,3,4) = / E Afc^Tr[r"-(i)r"-(2)r'^-(3)r"-(4)] [A^^^^{ail),a{2),ai3),a{A)) 

+ ^4mM1)'^(2)><^(3),^(4))" 
+ / E iVcTr[r'^-(i)r'^-(2)] Tr[r"-(3)r"-W] Ai.i^^{a{l), cj(2); cj(3), a{A)) , 

(2.4) 

where Ai-i^i and ^4;i,3 are 'partial amplitudes'. (Our fundamental representation color matrices 
are normalized by Tr[T"T^] = 5"''.) The notation '6*4/^4' denotes the set of all permutations of 
four objects 6*4, omitting the cyclic transformations. The notation '£'4/^2' refers again to the set 
of permutations of four objects, omitting those permutations which exchange labels within a single 
trace or exchange the two traces. That is, S4/Z2 = {(12 3 4), (13 24), (14 2 3)}. There is a similar 
decomposition for contributions linearly proportional to the number of particles in the fundamental 
representation, but with each term lacking one factor of Nc (see eq. (|3.1|) below). 
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These color decompositions have the virtue of making manifest the gauge invariance of the par- 
tial amplitudes. In addition the computation of the leading-color partial amplitudes from their cuts 
is straightforward. However, the cut structure of the subleading-color partial amplitudes is more 
complicated. For this reason it will prove convenient to use a somewhat different color decompo- 
sition, based on structure constants f"'^'^ and specialized to the all-plus helicity amplitude, whose 
specification we postpone until section ^ 

In evaluating the cuts we use the 't Hooft-Veltman (HV) variant of dimensional regulariza- 
tion [0] . (For a discussion of variants of dimensional regularization see refs. andrefs. 
Note however, that the one-loop near-equivalence of the four-dimensional helicity scheme and of di- 
mensional reduction assumed in refs. []T7| may not hold at two loops.) In this scheme observed 
external polarization vectors and momenta are four-dimensional, while unobserved and internal ones 
are D-dimensional. Cut lines are considered 'internal'. We use the convention that D > 4 (e < 0) so 
that a four-dimensional vector is given by setting the (— 2e)-dimensional components to zero. The 
HV variant violates supersymmetry Ward identities, but is commonly used and is straightforward 
to implement. 

In sewing the cut gluon lines we use the transverse, or physical state, projector, 

P..(^,r) = -r?,. + ^H!i_±^, (2.5) 

where £ is the gluon momentum and r is an arbitrary null 'reference' momentum that drops out of all 
final expressions. The Lorentz indices in this expression are all D-dimensional, but for convenience 
we will take r to lie in the four-dimensional subspace. For cases where only a single gluon crosses 
the cut, the second term in the projector drops out since it is contracted with conserved currents; 
for cases where two or more gluons cross the cut, the second term ensures that unphysical states do 
not circulate in the loops. We denote the number of gluon states circulating in the loop by Dg — 2, 
which is implemented by the Minkowski metric sum ri^'^r]^^, = Ds- In the HV scheme |14] which we 
are using here, Dg = D should be taken in the final expressions. However, it is useful to keep Dg 
and D independent in intermediate steps. 



3 Two-loop amplitudes with an internal scalar loop 

Although the double cut ( |2.3|) contains only tree amplitudes, the gluons crossing a generalized cut 
carry Lorentz indices along with them, making the algebra still non-trivial. Contributions with 
scalars crossing the cut are easier to compute first, and provide useful information about the gluon 
amplitude. Indeed at one loop, the contributions to identical-helicity all-gluon amplitudes with a 
scalar circulating in the loop are precisely equal, up to a trivial overall constant, to the contributions 
with either a fermion or gluon in the loop . This may be understood in terms of supersymmetry 
Ward identities An{+, +,...,+) vanishes to any loop order in a (supersymmetrically regulated) 
supersymmetric theory, and the combinations of gluon-|-fermion and scalar+fermion circulating in 
a single loop are supersymmetric. 

At two loops, even were we to use a manifestly supersymmetric regulator, we should not expect 
the simple one-loop proportionality between scalar and gluon loop contributions to persist. Nonethe- 



5 



less, it is plausible that contributions from different internal particles will be closely related. Indeed, 
for the all-plus helicity amplitude, we have found it useful to evaluate first the contribution with an 
internal scalar circulating around one of the loops. It serves as an excellent guide to an ansatz for 
the pure gluon case. 




Figure 4: Representative diagrams for the contributions with a single scalar loop. The class (a) diagrams 
can be drawn on the plane with all external legs on the outside, and one internal gluon line on the inside of 
the scalar loop. The class (b) diagrams have all the gluon lines on the outside of the scalar loop. 



To illustrate the utility of scalar-loop contributions in constructing the gluon-loop one, consider 
the two- loop four-gluon amplitudes with a single scalar loop. Figure |^ shows that there are two 
different classes of planar diagrams, divided according to whether or not an internal gluon line runs 
inside the scalar loop. For the case of scalars in the fundamental representation of SU{Nc), the 
contributions proportional to the color trace Tr [T^^T^^T^^T""'] are 



Tr [r'^i 2^^2 2^<i3 q^ai j 



2- loop, (b) 



1 

'Nr. 



A 



2-loop,(a) 



(3.1) 



where Nc is the number of colors, and the labels (a) and (b) correspond to the diagram classes 
depicted in fig. ^. The amplitudes are color-stripped; that is, they should be computed using color- 
ordered Feynman rules from which color factors have been removed. Because the diagrams of 
class (a) and (b) enter eq. ( |3.lD with different powers of Nc, they are separately gauge invariant and 
may be computed independently. For scalars in the adjoint representation the two contributions are 
still separately gauge invariant, but there is no relative factor of —N^- The leading-color contribution 
proportional to the same color trace is then 



^2-loop,(b) ^ ^2-loop,(a) 



(3.2) 



We begin our evaluation of scalar-loop contributions with the s-channel double two-particle cut 
for the subleading-color term in fig. §(a). This is given by the product (|2.3| ) where only scalars cross 
the generalized cuts, 



^2-loop,(a)^^+^2+,3+,4+) 



2x2- cut 



f-2) 



2 \2 



[.4; ^ ^^^°'^{il,i^,3~^ ,4^^) 



X A 



[12] [34] / 



•S12A1A4 



I 

+ 7:- 



Af A4 



(3.3) 



where Sij = {hi + kj) . The vectors Aj represent the (— 2e)-dimensional components of the loop 
momenta ii; that is, ii = + Aj, where 4^' has only four-dimensional components. Superscripts ± 
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label helicities of external gluons, while a superscript s indicates that a cut leg is a scalar. (A cut 
momentum without a superscript will implicitly denote a gluon.) We used the color-ordered Q tree 
amplitudes (or partial amplitudes), 

<^*^(l+,2+,-£i,-£J) = iA? ^^^^ 1 



Atieeffs ps _ps _fs\ 



(12)(4-fci)2' 
[h - 2 



^4 ''"(^4, ^3,3+, 4+) can be obtained by relabeling ^*/^^(l+, 2+, -^i, -£f ). We have expressed the 
amplitudes in terms of spinor inner products [l| which are denoted by (ij) = and [ij] = 

where are four-dimensional massless Weyl spinors of momentum k^, labeled with the 
sign of the helicity. 

As a working hypothesis, from eq. (^]^) we take the subleading-color planar amplitude to be[| 



(3.5) 



.2-loop,(a) „+ + ^ — Tc TPr\2\2i/ „ ^ _i_ ''^ rbow-tic r a, 2 \ 2 1 1' „ ^ 

^4 (i ,2 ,3 ,4 ) - l-^J^y^^[Si2l4[Ap\\[Si2,S23) + -I4 [Ap\\[Si2) 

+ (Sl2 ^ S23) 

where the two terms correspond to the two terms in the second equation of eq. ( |3.3| ), 

'^^['P{>^i,P,q,ki)]{si2,S23) 

d^'p d^q V{Xi,p,q,h) (3.6) 



(27r)^ (27r)^ p2 ^2 + ^)2(p _ ^^)2 (p - - yta)^ - ^4)^ {q - ks - k^f 

is the planar double box integral displayed in fig. ^(a), and the vectors Ap, \q represent the (— 2e)- 
dimensional components of the loop momenta p and q. The numerator factor V{Xi,p,q,ki) is a 
polynomial in the momenta. The 'bow-tie' integral X^"™"*"^ shown in fig. ^ is defined by 

l'^"'^-'"'[V{K,p,q,ki)]{si2) 

d^^p d^q V{K,p,q,ki) (3.7) 



{2tt)D (27r)^ p2 ^2 _ ^^)2 (p-ki- k2f [q - k^f {q-k^- k 
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\2 



The one 'missing' propagator, compared with if , makes X4°^"*"^ much simpler to evaluate, since it 
is a product of two independent one-loop integrals. Although it is not manifest, the overall spinor 
prefactor is symmetric under all 4! permutations of the external legs, including the cyclic relabeling, 
^1 — > ^2 — *■ ^3 ^ ^4 — > ^1. Hence the amplitude (^), including the term with S12 and S23 
exchanged, is cyclicly invariant too. The exchange term generates the proper double two-particle 
cut in the t-channel. 

Cutting a line forces the corresponding propagator to be present, uncanceled, in all integral 
functions. In general, an ansatz for an amplitude based on a limited set of cuts will fail to capture 
integral functions which are missing some or all of the cut propagators. As one refines the ansatz by 
investigating additional cuts, it may be necessary to add in integral functions that do not have cuts 



^These bare amplitudes should be multiplied by the scale factor {fJ-o)^'^, and after renormalization instead by the 
scale (Hr)^'^- Here we take these scales to equal unity except where noted. 
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1/ \4 

(a) (b) 

Figure 5: The planar and non-planar double box integrals. 




Figure 6: The 'bow-tie' integral appearing in eq. (3.5) 



in previously-investigated channels. For example, the integral depicted in fig. |^ contains only three- 
particle cuts and would not be detected in two-particle cuts. However, for the planar amplitude 
under consideration, it turns out that such integrals do not appear, and the double two-particle cuts 
are in fact sensitive to all functions appearing in the amplitude. 




Figure 7: An integral function containing only three particle cuts. 



To verify that eq. (^]^) is indeed the complete amplitude corresponding to the class (a) diagrams 
in fig. ^ we must evaluate the two- and three-particle cuts. To evaluate the two-particle cuts we 
need the first of the tree amplitudes given in eq. ( |3.4D , as well as the subleading-in-color one- loop 
amplitude (see fig. ^(a) for a representative diagram). 



^l-loop,(a)^^+^2+,3^4^) = -M 



(12) 



512 2", 



1-loop 



1 



Ap](si2,S23) + ^2:3 



l-l00p,(4)r^ 21/ N 



(3.8) 



(We compute all required one-loop amplitudes from their s- and t-channel cuts.) We have taken the 
momenta of the external legs 1 and 2 to be four-dimensional while the momenta of legs 3 and 4 are 
D-dimensional, so that they can be inserted into the two-particle cuts. The one-loop box integral in 
eq. (^) is defined by 



^l-loop 



[P{Xp,p,ki 



V{Xp,p,ki) 



(27r)^ p2 _ ^^)2 (-p _ /.^ _ ^2)2 (p + kif 



(3.9) 
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Figure 8: Representative planar diagrams for the subleading-in-color (a) and leading-in-color (b) one-loop 
amplitudes for two gluons and two scalars. The diagrams are divided into the two classes according to whether 
or not there is an internal gluon line to the right of the scalar line. 



The triangle integral Tg ^""P'^"^-* is obtained by removing the l/{p + k^)"^ propagator from eq. (jS 
Taking the two-particle cut of the ansatz (^]^) we obtain 



A 



2- loop, (a) 



(l+,2+,3+,4+) 



^-^°°P'(")(^|,ri,3+,4+) 



2-cut ^ . ^ ^ ^. , ....... 

^^i-ioop,(a)^^+^2+,-£^,-£^) X AT^£iel,3+,4+), 

as required. The two-particle cuts work simply in this example because all terms in the one-loop 
amplitude ( |3.8D could be detected via the s-channel cut. 

The three-particle cuts are more intricate, but may be handled numerically. For the s-channel 
case, we need the product of tree amplitudes. 



(3.11) 



where fi, v are Lorentz indices for the single gluon line crossing the cut, and the physical state 
projector P^^y is defined in eq. ( p.sp . If our ansatz ( |3.5| ) is correct, then the cut expression ( p.ll ) 
should be equal to 



. [12] [34] 
'(12) (34) 



Ag 



S23 



+ 
+ 



Sl2 



{h + ^2)2(^3 - fe2)2(^2 + ^3)2(^1 + k^)^ 
Sl2 



(3.12) 



{h + £2f{h + ^?,)\t2 + ^3)2(^1 - kxf_ 

which is obtained by cutting the two planar double box integrals in eq. (^]^) as shown in fig. |^. The 
bow-tie integrals do not contribute to three-particle cuts. The equality of these expressions should 
hold before any loop or phase-space integration is performed. 

The expressions for the cuts are smooth analytic functions with no explicit dependence on the di- 
mension D of the loop-momentum vectors. It is therefore sufficient to verify the cuts for integer values 
of -D, with > 4, where numerical evaluation is straightforward. We randomly generate D-vectors 
satisfying the appropriate cut kinematics — in this case, seven null-vectors A;2, /ca, A:4, ^i, £2, ^3} 
obeying 

kx + k2-lx-l2-h = k3 + k4 + £i+h + h = 0, (3.13) 

where the ki are four-dimensional, i.e., they have vanishing extra-dimensional components. We then 
check whether the cuts agree with our ansatz at this point in phase space, in this case by comparing 
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1 



1 



Figure 9: The s-channel three-particle c uts o f the class (a) scalar contribution in eq. (3^), and of the 
leading-color pure gluon contribution in cq. (4.1). The dashed lines represent the cuts. 



numerical values for eqs. (|ll|) and §J^). Because no integration is required, the check can be 
performed to very high accuracy (25 or more digits), and repeated quickly for other phase-space 
points. As mentioned previously, we take the dimension Dg, which appears in the number of gluon 
states circulating in the loop, to be independent from the dimension D of the loop momentum 
vectors. To distinguish between e.g. (Ai • A2)^ and Af A^ (the most subtle degeneracy arising at two 
loops) it is necessary to choose D > 6. 

Together with the two-particle cuts, the numerical verification of the three-particle cuts proves 
that the ansatz ( p. 51 ) does indeed give the complete subleading-color contribution to the color coef- 
ficient Tr[r"ir''2 7-"3 j^""*] for the case of a scalar in the loop. 

Following a similar procedure we have also obtained the results for the leading-color contributions, 
obtained from the diagrams of the type displayed in fig. |^(b). The color-stripped amplitude for case 
(b) is. 



^2-ioop,(b)(^+ 2+,3+,4+)-^-ii^ 



suI^iXl Ap+„](S12, S23) 



(12) (3 4) 

+ ^^^^T4"[A^ X'^ Hp + qf + .12)] (.12, .23) (3-14) 



+ 4T4^°"-t'^[A^ (Ap • A,)](si2) + cyclic perms of (1, 2, 3, 4) 



S12 

rbow-tie r \ 2 



In this case the number of gluonic states {Ds — 2) circulating in the loop appears. Although some 
of the terms in eq. ( 3.14| ) vanish under integration, we have kept them because they do not vanish 



in the cuts before integration and are therefore important for verifying the cuts numerically. 



(b) 

Figure 10: Representative planar diagrams with two scalar loops. The ones in (a) have two distinct loops 
connected by a gluon line; the ones in (b) have a quartic scalar interaction. 



The two-loop four-gluon amplitude with scalars present can also receive contributions with two 
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scalar loops. Figure [l^(a) shows a class that is present just due to the gauge interaction. Because 
supersymmetric gauge theories with scalars always include quartic scalar interactions (arising from 
I?-terms), it is useful to also calculate the class shown in fig. lC(b). Both types of contributions are 
simple, 

• double— scalar, (a) -|- q+ a + \ . [12] [34] 1 

Tbow-tier\2 \2 1 1 ^ , „\2 , 1„ M/„ „ \ 
^4 (J- ,^ ,-3 ,4 j OC I ^^^^ ~2 [\\v(P + 9j + 2^12jJ(,Sl2,S23j , 



^4 '^"^(l+,2+,3+,4+) oc iJ^^I|^Tfw-uc[^.^.](^^^)_ 

(3.15) 

The precise normalization is unimportant; in constructing the pure gluon ansatz we shall add these 
terms to eqs. ( t3.5| ) and ( |3.14| ) with coefficients to be determined later. 



4 The pure gluon amplitude 

We now turn to the two-loop pure gluon amplitude ^4 '°°P(1+j 2+, 3+, 4+), which is the main subject 
of this paper. It is convenient to begin with the leading-color contribution generated by planar 
graphs. 

4.1 Leading color 

The planar amplitudes with internal scalars from section ^ serve as a guide toward constructing the 
planar pure gluon contributions. We write an ansatz in terms of loop integrals by taking a linear 
combination of the above results for one scalar in the loop, plus contributions with two independent 
scalar loops. This turns out not to be the complete answer so an additional polynomial in the Aj 
is required, whose form is not too hard to guess. We numerically evaluate the Z)-dimensional cuts 
at a number of different phase-space points, and use this information to solve for the unknown (but 
kinematically constant) coefficients of possible additional terms. In this way we obtain a simple 
representation of the planar pure gluon amplitude in terms of loop integrals. The form of the planar 
amplitude then suggests an ansatz for remaining non-planar contributions, which we again verify 
numerically from the cuts. 

In somewhat more detail, our starting point is the contribution of a single adjoint scalar loop, 
eq. (|3.2| ). We multiply this contribution by the number of gluon states, Ds — 2. This would produce 
the complete answer at one loop, but here we must make two corrections: 

• a factor of 1 /2 for terms proportional to {Ds — 2)^ , from different combinatorics of contributions 



in which the two loops have no common propagator, as in fig. 10(a). 



• the addition of a term (with an unknown coefficient) induced by a quartic scalar interaction, 
as in fig. pX|(b). 

Even with these corrections, the ansatz does not quite work. However, it does satisfy the double 
two-particle and three-particle cuts in D = 4 (i.e., when all (— 2e)-dimensional components are set 
to zero) and also in D = 5. This provides a strong clue to the form of the additional term; it must 
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be proportional to powers of Aj, and it must vanish when all the vectors Aj are parallel. The missing 
term, determined numerically, is the double box planar integral with numerator 16{Ds — 2)si2(^{Xp ■ 

Ag)^ — Ap Aq^ . We then find for the pure gluon leading-color all-plus helicity amplitude in eq. (|2.4|), 

4?^i(l+,2+,3+,4+) =^^34 + ^2341, (4.1) 

where the 'primitive' planar amplitude ^1234 is analogous to the one appearing in the = 4 super- 
Yang-Mills case 0, 



^1234 



[12] [34] 
(12) (3 4) 



•7-P 



\Ds - 2){\l \l + \l + \l + 16((Ap • \f - \l A2)] {Si2^ S23 



+ A{Ds-2)Z 



bow-tie r 



+ 



S12 



T-bow-tie 



XlXliip + qf + S12) isi2,S23)}. 



(4.2) 



4.2 Subleading color 

We may observe that all terms in A^*^ have two-particle cuts in exactly one channel. If we assume 
that this property also holds for the subleading-color contribution, then performing the sewings 
implied by eq. (|2.l| ), we find that full two- loop amplitude can be written in the following form 



^f'°°P(l+,2+,3+,4+) 



C1234 ^i'234 + C'3421 ^3421 +^"12^34 ^12^34 + ^'34^21 ^34^21 + C(234) , (4.3) 



NP 



where '+ C(234)' instructs one to add the two non-trivial cyclic permutations of (2,3,4). The values 
of the color coefficients may be read off directly from fig. ^. For example, Cf234 is the color factor 
obtained from diagram (a) by dressing each vertex with an f"'^'^, where 



abc 



iV2f 



abc 



Tr([r'',r*]T'= 



(4.4) 



and dressing each internal line with a 6"''' . Similarly, C^-^^^ is obtained by dressing diagram (b) . The 
symmetries of the color factors can also be read off the diagrams: 



C 



a 



1234 — "-"3412 - 

/^NP _ ^NP 
*-"12;34 — '-"21;34 



a 



2143 
^NP 



a 



4321 5 
-yNP 



'-"12;43 — '-"21;43 • 



(4.5) 



The corresponding planar and non-planar primitive amplitudes share the same symmetries with their 
associated color factors. Due to these symmetries, ^°"P(^i+ ^ 2"*", 3+, 4+) has the required total (S'4) 



permutation symmetry, even though only six permutations appear explicitly in eq. ( |4.3| ). 

In the amplitude ( |4.3| ) no ultraviolet subtraction has been performed. The subtraction in the 
MS scheme is given below in eq. ( ^.21 ). 

The one-channel assumption about the two-particle cuts is non-trivial and is not expected to hold 
for general helicity amplitudes. This property and the decomposition of eq. ( |4.3| ) are also satisfied 
by the A'^ = 4 supersymmetric amplitude ]|^, ^. One may verify the validity of the assumption a 
posteriori by checking the two- and three-particle cuts, eqs. (p.ip and (|2.2[). Although it might appear 
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that the expression ( |4.3D does not assign the proper color factors to bow-tie integral contributions, 
the unwanted terms (the subleading-color parts of C^) cancel in the permutation sum. The partial 
amplitudes appearing in the standard color decomposition ( |2.4D can be expressed in terms of ^4^ and 
A^P. The equations are identical to those for the = 4 supersymmetric case discussed in ref. m. 

Prom the point of view of the two-particle cuts, the planar amplitude arises from sewing a four- 
point color-ordered tree to a color-ordered one-loop amplitude, where the sewn legs are adjacent in 
both amplitudes. In other terms in eq. the permutation of legs on the one-loop amplitude 

A\ results in sewing non-adjacent (diagonally opposite) legs to the tree. This gives rise to a 
primitive amplitude built out of non-planar two- loop integrals, 

^12^34 = I 2 1 Li si^lf^ i^s - 2)(Ap + Xp Ap+g + Xg Ap+g) + 16((Ap • A^)^ - A^ A^) (si2, S23) , 



(12) (34) 

where the non-planar double box integral, depicted in fig. ||(b), is given by 
lf^[V{Xi,p, q, ki)]{si2, S23) 



(4.6) 



d^p d^q V{Xi,p,q,ki) (4.7) 



{2tt)D (27r)« p2 g2 + g)2 _ fc^)2 _ ^2)2 {p + q + k^y [p + q + k^ + k^f ' 

Note that ^12^34 is symmetric under ki k2, and independently under k^ /C4. The non-planar 
amplitude has no bow-tie contributions; however, the term proportional to ^(Ap • A^)^ — A^ A^^ 
persists. 

To verify that the expressions in eqs. (^]^), ( |4.2D and ( |4.6| ) give the correct full gluon ampli- 
tude, we have evaluated all the independent cuts. Consider, for example, the leading-color planar 
contribution proportional to the color factor iV^ Tr[T"iT"2T"3T"''], which is given by eq. (4.1). For 
this contribution, the si2-channel three-particle cut integrand should be equal to the product of 
color-ordered five-gluon tree-level partial amplitudes, 

4-^(1+, 2+,-£g, -0 xP^,(^i,r)P,^(^2,r)Pp-,(£3,r)x4'--(^?,£f,£2,3+,4+). (4.8) 



We have verified numerically that the product of tree amplitudes (4.8) is equal to the three-particle 
cut of the leading-color expression ([4.1|), given by 

..[12] [3 4] 



(12) (34) 



{Ds - 2){XlXl + XlXl + AlAi) + 16((Ai • Aa)^ - AfAf 
fl2 , . 

kiY{i3-k2Y{h + k3Y{h + k^Y^ {i^ + i2Y{h-k2Y{i2 + h?{h + k^Y ^ • ^ 

S12 



+ 



(^1+ ^2)2(^3 + ^3)2(^2+ 4)2(^1 -fei: 



The three different terms in the last factor in eq. (|4.9| ) originate from the two cut planar double box 
integrals in the same way as in the scalar example ( |3.12| ), as illustrated in fig. 

As a second example of a si2-channel three-particle cut, but one that is sensitive to non-planar 
(subleading-color) contributions, we evaluate the product 

4'--(l+,-£g,-£^,2+,-<) X P,c.{h,r)P,p{t2,r)Pp,{h,r) x 4'--(^?, £f , 3+, 4+) (4.10) 
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numerically and find that it equals 



..[12] [34] 



(12) (34) 



{Ds - 2){XlXl + AfAi + A^A^) + 16((Ai • As)^ - A^A; 



Sl3 Si3 



ih - h)Hh - k2)^l3 + k3)Hei + {£3 - kif[l2 - k2Y{l3 + k3?(h + kif 
_^ £23 ^ S12 



{h - kiY{i2 - k2Y{l3 + k3)Hh + k^y {£3 - kiY{i2 - k2Y{£2 + i3Y{h + k4 
S12 S12 



(ii - kif{l2 - k2f{l3 + k3Y{h + I2Y ih - kiY{£2 + + k3?{li + £2? 

S12 



{h-k2Y{l2+£3?{£3 + k3?{£l+£2)\ 

(4.11) 

The latter candidate expression is found by extracting from eq. ( [4 .31 ) the seven terms with the 
correct color ordering to contribute to the cut ( [4.1[1| ). (Three different planar double box integrals 
and three different non-planar ones appear; one of the non-planar integrals can be cut two separate 
ways. Minus signs are due to the antisymmetry of the structure constants in the color factors.) 
More generally, we have verified that the amplitude (4.3) has the correct D-dimensional two- and 
three-particle cuts in all channels, proving it to be correct. 

We have also shown that the double two-particle and three-particle cuts of the Feynman gauge 
Feynman diagrams, including the ghosts, match the cuts of our expressions. In performing this 
consistency check there is no need to include physical state projectors on the intermediate gluon lines; 
the ghosts automatically cancel the unwanted longitudinal modes. We carried out the comparison 
numerically for specific values of the (cut) loop momenta. This provides a non-trivial check that 
our expressions are in one-to-one correspondence with the results that one would obtain via a direct 
evaluation of the Feynman diagrams. 

In the appendix we perform the final step in a closed-form evaluation of the two-loop amplitude: 
We present the values of the integrals appearing in the amplitudes, expanded in e through 0{e^) 
and expressed in terms of polylogarithms. The bow-tie integrals are rather simple to evaluate since 
they are just products of one-loop integrals. The planar and non-planar double box integrals are 
much more difficult to obtain, because they have poles up to order 1/e^ and an intricate analytic 
structure (especially the non-planar cases). 



5 Divergences 

In this section we compare the divergence structure of the all-plus helicity amplitude ( |4.3| ) against 
the one expected from general principles. Catani has previously presented a universal factorization 
formula for the infrared divergent parts of dimensionally regulated, renormalized two-loop ampli- 
tudes [12|. In the all-plus case, the divergence structure is relatively simple, due to the vanishing 
of the corresponding tree-level helicity amplitude; both the ultraviolet and infrared divergences are 



essentially the same ones encountered at one loop [13|. Nevertheless, we find it convenient to use 
Catani's formalism, because of its more general applicability. 

In the color space operator language used by Catani, the infrared divergences of a renormalized 
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two-loop amplitude are, 

+ /i(6,/x|;M) (^|;M))H.a + finite, 



(5.1) 



where |7Wi^^(/i|; M)) Rs. is a color space vector representing the renormalized L loop amplitude. 
The subscript rs stands for the choice of renormalization scheme, and /Ur is the renormalization scale 
(which we have set to unity elsewhere in the paper). These color space vectors give the amplitudes 
via, 

A(l''\---,n"") = {ai,...,an\Mn{pi,...,Pn)) , (5.2) 

where the Oj are color indices. The divergences of An are encoded in the color operators I^^^ (e, /x^; {p}) 
and I^'^\e, /x^; {p}). For the all-plus helicity case, the tree amplitude vanishes so I^'^\e, {p}) does 
not enter. The operator I^^\e, fij^; {p}) is given by 



2 {An 



2-e 



Til 



i=l 



ri /;u|e-^--y 




.e^V 2pi-pj J 


e\ 



(5.3) 



where Xij = +1 if i and j are both incoming or outgoing partons and Xij = otherwise. The 
color charge Tj = {T/*} is a vector with respect to the generator label a, and an SU{Nc) matrix 
with respect to the color indices of parton i. For the adjoint representation = if^"-^, and so 
T? = Ca = Nc- Also, for external gluons in a pure glue theory, = ^Ca- (We have not included 
an overall factor of e~'^^*^^^ appearing in Catani's expression, but have included a factor of 2/(47r)^~'^ 
because of differing overall normalization conventions.) 

Specializing I^^\e, fij^; {p}) to the pure gluon case, and then rewriting the divergent terms in 
eq. (^) in our notation gives. 



^2-loop, ren.|--|^ai ^an' 



div. 



^^(ij),2-loop,ren.(-^^2,3,..., 

i<j 



n, 



div. 



where 

__4(i,i),2-loop, ren. ^-^^ 2, 3, 



div. 



-g erf ' ' f 



biaiC ?cajbj 



1 



+ 



11 1 

:b- 



and f is defined in eq. M 



X ^^-'°°p(r\2''2, 

In converting the overall normalization we used 
1 r(l + e)r2(l-e) 



cr 



(47r)2-^ r(l - 2e) 



(47r)2-^ r(l - e) 



(5.4) 



(5.5) 



(5.6) 



Equation ( ^ ) holds for every pure gluon helicity amplitude of the form ^„(ib, +,+,..., +); 
such amplitudes vanish at tree-level due to supersymmetry Ward identities ||ll| . (Although the pure 
gluon theory is not supersymmetric, its tree amplitudes obey supersymmetry identities because the 



15 



gluons' fermionic superpartners cannot contribute at tree level.) Specializing now to the four-point 
case, we have 



_^2-loop,rcn.^^^2,3,4)_^. ^ ^ ^, 
div. — ' 



I (j j),2-loop, ren. 



(1,2,3,4) 



i<j 



div. 



(5.7) 



We shall evaluate the term in eq. (5.7) with = (1,2), 



div. 



I, 111 



jbiaic jca2b2 



6 e. 

^^"'°°p(l^^2^^3''^4»*), 



(5. 



and then obtain the other five cases by relabeling i and j. Note that we count (i, j) = (3, 4) as being 
distinct from = (1,2) in the sum { f).7\ ). 

The one-loop amplitudes may be decomposed in terms of SU{Nc) structure constants ||2^ (see 
also eq. (2.30) of ref. [|), 



4-^°°P(l, 2, 3, 4) = ^r°'^(l, 2, 3, 4) + ClZ"" <"°'(1, 2, 4, 3) 



yl-loop A 1-loOp 



1-loop A 1-loop / 



+ C^-^°3°P4-^°°P(1,4,2,3) 



(5.9) 



where Cj^234''' color factor obtained by dressing each vertex in the one-loop box diagram in 



fig. 11 with a structure constant / , and dressing each bond between vertices with a 6 . In 



this decomposition the A\ are free of all group theory factors, and are invariant under cyclic 
permutations of their four arguments. 

For the all-plus helicity case, the primitive amplitudes A] are proportional to a one-loop box 



integral [plj . 



yir^°°P(l+,2+,3+,4+) = {D, - 2) MM T]-i-P[A^](.i2,.23) , 



with 



4"'°°'^[A^](^12,S23)^ / 



(2^ p^{p-kif{p -ki- k2?{p + ki) 



2 ■ 



(5.10) 
(5.11) 



Figure 11: C'^234°^ given by dressing each vertex in the figure with /"''^, and dressing each bond between 
vertices with a 6°^''. 



Using eq. (|5.9|) the divergence (|5.8|) is 



A 



(1,2), 2-loop,ren. 



1,2,3,4) 



div. 



I, 111 

6 e 



Jb,a,cjca,b, (^^1-loop ^l-loop^^^ 3, 3, 4) 



+ CltZ 2, 4, 3) + C™ Ar^°°P(l, 4, 2, 3) j . 

(5.12) 



1-loop < 1-loop^ 



16 



X 



-3 b2 

+ 

-4 h 



'4 4- 

+ 

~3 h 



Figure 12: A schematic version of the color factors in eq. ( p. 15 ). The color factors are given by dressing 
each vertex with Z"''^, and dressing each bond between vertices with a 5"''. 



We display a schematic version of the above color factors in fig. 12. Anticipating our comparison 



to the divergences appearing in the two-loop amplitude (O), we observe that the one-loop color 
coefficients can be related to the two-loop ones via, 



biaic /ca2fe2 /^t"l°°P 



joiaic J, 



6lfe2 34 



c 



/ biaic f £02^2 /^t-loop 



J J '-"6i4b23 ~ '-"34;12 ■ 



Using these relations in eq. (|5.8| ) we have 



1234 ) 

■yP 

'12' 

-yNP 



(5.13) 



^(l,2),2-loop,ren.^^^2,3,4) 



div. 



1 



11 1 
¥7 



(^^f234 4"'°°''(l,2,3,4) 



+ Cf,,s ^^"''(l, 2, 4, 3) + CZ2 ^r^°°'(l, 4, 2, 3) 



(5.14) 



so the color factors are the same ones appearing in the two- loop amplitude. Inserting eq. ( 5.14| ) into 
eq. (^), summing over i and j, and making use of eq. (|4.5|), we obtain the predicted divergence, 



^2-loop,ren.(-^^2,3,4) 



div. 
-S12) 



+ 



11 1 



6 e 



1 l-loop / 



(2CP 344-^°°P(l, 2, 3, 4) + 2C3^42i^r'°°'(3, 4, 2, 1 

+ 4-^°°P(l, 3, 2, 4) + CZi Al-'°°^(3, 2, 4, 1)) + C(234) . 

(5.15) 

Now we compare eq. ( ^.15 ) to the divergences found in the two-loop all-plus helicity amplitude 
( |4.3| ), using results for the two- loop integrals given in the appendix. All of the integrals are finite 
except for T4'[Ap_|_gAg], 2^4'[Ap_,_gAp] and ^^^^[ApAg]. For eq. ( |5.15| ) to be recovered, it is crucial that 
the divergent parts of these integrals are proportional to the same one-loop box integral ( ^.11 ) 
encountered in A] 



^4 [■^p+g'^g]('5l2,S23; 



•p+q^q 



2^4 [-^■p+(;-^p]('5l2,S23^ 



div. 



J] 



rNP r \ 2 x 2 



div. 



div. 



-icr (-S12)' 



1-e -^l-loop |-\4 



[AJ(S12,S23) , 



-«cr^(-si2) ^ 'x]"'°°'^[Ap](si3, 



S23j 



(5.16) 

as can be seen from eqs. ( |A.ll ) and ( A.14| ). 

This factorization of the singular behavior of the integrals can be understood heuristically. The 
1/e^ poles come from loop momenta that are simultaneously soft and collinear with two adjacent 
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external legs, so that three consecutive propagators can go on shell. Because the external momenta 
are strictly four-dimensional, the singular region of loop momenta is also four-dimensional. However, 
a factor of Xp in the numerator suppresses the region where p becomes four-dimensional. Thus in 
X|'[Ap^gAg], for example, the region with ~ and ~ dominates the integral, while q and 
p + q remain off shell. The space-time picture of the dominant integration region corresponds to 
shrinking the rightmost box of the planar double box in fig. ^(a) to a small size. The kinematics of 
that inner box integral is hardly affected by the soft momentum p, so it is essentially the same as 
the external kinematics, and it can be factored out of the outer triangle integral that remains. The 
triangle integral gives the factor of — icr (— •si2)~^~'^/e^ 'dressing' the box integral. 

Apart from those in eq. ( |5.16| ), all the other integrals appearing in the two- loop all-plus amplitude 
are finite in the limit e — > 0. For these integrals, whenever there is a loop momentum from the set 
{p, q, p + q} whose four-dimensional limit is not suppressed by a numerator factor of Xf, then that 
loop momentum has at most two propagators containing it, not three. This property prevents an 
infrared divergence from occurring. 

Using eqs. ( p^ ), ( |42D and (|^, and the total symmetry of the prefactor [12] [34] /((I 2) (3 4)), 
the divergences of the bare (unrenormalized) primitive amplitudes are 



^1234 



^12:34 



= -2 cr ^(-.i2)-^4''°°'(l+, 2+, 3+, 4+) , 
div. e 

=-cr^(-si2)-^4''°°'(l"',3+,2+,4+). 
div. e 



(5.17) 



Inserting these expressions into the bare full amplitude (4.3) then yields 
^f^°°P(l,2,3,4) 



^. = -g cr 

div. 



|(-si2)-^Cf234^1-'°°P(l,2,3,4) + |(-si2)-^C3V^^-^°°P(3,4,2,l) 

+ ^{-su)-'Cf,% ^ri°°P(l, 3, 2, 4) + li-sur'Cl% 4-'°°P(3, 2, 4, 1) 
-f C(234) . 

(5.18) 

The expression (|1|) already agrees with eq. ( |5.15| ) at the level of the {—Sij) ^/e^ poles. To 



check the agreement at the 1/e level, we must first renormalize our two-loop amplitude in the MS 
scheme. The MS counterterm to be subtracted from '°°^(1, 2, 3, 4) is 

C.T. = 4/cr^-4-^°°P(l+,2+,3+,4+), 

= 4A^,/cr^^[c;2t4°"4''°°"(l,2,3,4) +C3\tr • 

(5.19) 

The relative simplicity of the ultraviolet subtraction term is due to the vanishing of the corresponding 
tree-level helicity amplitude; the counterterm has the same structure encountered in other helicity 
configurations at one loop, up to the overall factor and the replacement of a tree amplitude with the 
one-loop amplitude. (At two loops, an additional 1/e^ subtraction is required for helicity configura- 
tions with non- vanishing tree amplitudes. Additional subtraction terms are necessary for amplitudes 
including fermions or scalars as well.) 
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Finally, we use the color conservation identity J^i'^i = l^^l) which implies that 



n 



Nc \Mn) = E 1-^") = - E • 1-^") ' 



(5.20) 



when all external particles are in the adjoint representation. Replacing the ANc in eq. ( ^.191 ) with 
the help of eq. ( 5.20| ), the color structure of the MS counterterm (and its dependence on A\ 



can be put in exactly the same form as Catani's expression for the divergent terms, eqs. (5J-), (|5.3| ) 
and (^.4|). We can then process it into the form ( |5.15| ) in just the same way, obtaining 



C.T. =/cr 



11 1 
( 

3 e 



^ 344"'°°'^(1, 2, 3, 4) + y ici; 2i4-i°°P(3, 4, 2, 1) 



+ ^ -^CZ, 3, 2,4) + ^ -^CZi ^4"'°°''(3, 2, 4, 1) 



(5.21) 



+ C(234) 



After including the ultraviolet subtraction ( p.21| ), the infrared divergences of our renormalized am- 
plitude agree perfectly with the expected form ( [5.15| ). 



6 Conclusions 



In this paper we presented the pure glue contribution to the identical helicity two-loop four-gluon 
amplitude in QCD. We calculated the amplitude using its analytic properties in D dimensions. More 
generally, any amplitude in a massless theory can in principle be obtained in this way. Although 
QCD is not supersymmetric, we also used supersymmetry as a helpful guide in constructing the pure 
gluon amplitude from amplitudes with scalar loops. Cases with quarks will be discussed elsewhere. 

We expect the type of methods used here to be useful for obtaining more general two-loop am- 
plitudes relevant for high energy processes. Although much more remains to be done, we are hopeful 
that calculations of next-to-next-to-leading-order multi-jet cross sections will become possible. 
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A Integrals 

In this appendix we evaluate the integrals appearing in the all-plus helicity amplitude, in an expan- 
sion in e through 0{e^). 

The planar bow-tie integrals are rather simple to evaluate because they are products of one-loop 
integrals, 

T4"[A^A^](.) = -i^, 
X^-*-[A2A2(p + g)2](.,t) = __L_i_(t -4.), 
T4'°"-*'^[A^(Ap-A,)](s)=0. 
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Two of the planar double box integrals and one of the non-planar ones vanish through 0{e^), 

l![XlXl]{s,t) = 0, 

ir[iXp-x,f-xlxl]is,t) = o. 

The non- vanishing planar double box integral, evaluated in its Euclidean region where s,t < 
and n > 0, is, 



^I[Xl+qXl]is,t) — j|'[Ap+gAp](s,t) 



C2r (-s) 



+ 



X 



1 1 

■- + - 

e 



^3 

where x = ^/s and 



■ Li3(-X) + C(3) + Inx Li2(-x) + 



-ln(l + x) ~ 3 ^"^^ ~ 2^"^ 3 



-IT 

~6 



Inx 



6 1 + x 



52 



The poly logarithms are |22] 
Li2(x) = - 



C2T 



dt 



1 r(l + 2e)r3(l-e) 



- (47r)4-2e 

ln(l - t) 



r(l - 36) 



t 



Li3(x) 



dt 



Li2(t) 



(In^x + TT^) 

(A.2) 
(A.3) 

(A.4) 



and (^(3) = Li3(l) = 1.202.... Equation (A.2) was originally found by direct integration over 
Feynman parameters, after performing certain subtractions to allow expansions of integrands in e. 
We have also checked that the same result can be obtained via the general tensor integral reduction 
method of Smirnov and Veretin . 

The remaining planar integrals in the full amplitude ( |4.3| ) may be obtained from eq. ( [A.2| ) by 
relabeling the legs, which induces transformations on X; for example. 



(A.5) 







X ^ 


1 

X 


t <-> li 




X ^ 


-1-x 


s ^ u 




X ^ 


-X 

1 + 



using s + t + u = 0. 

For two reasons, it is useful to represent the poles in e of eq. ( [A.2 ), and of the non-planar integral 
^^^[XpXq], in terms of the one-loop box integral ll ^°°^[Ap] defined in eq. ( ^.11| ). First, the divergence 
structure of the two-loop all-plus helicity amplitude becomes more transparent, as discussed in 
section ^ given that t1"^°°P[A^] appears in the corresponding one-loop amplitude. Second, the 
description of the finite {0{e^)) remainder terms is simplified, especially regarding their values in 
different kinematic regions. The heuristic reason why the singular parts of the two- loop integrals 
are proportional to X4 [A^] was given in section ^. 
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The one-loop box integral ( ^.11 ) can be performed by first integrating over the (— 2e)-dimensional 
components of p, which results in an integral of the scalar type but in 8 — 2e dimensions, not 4 — 2e. 
A standard one-loop 'dimension-shifting' identity allows one to trade this integral for a quantity 
C{s,t) which is essentially the scalar integral in 6 — 2e dimensions [^3|. This integral has neither 
ultraviolet nor infrared divergences, which makes it straightforward to expand in e. 

Concretely, we have 



rl-loOp r -i 4 



A^](s,t) = -6(l-e)T, 

-6(1 



1-loop, i:)=8-2ei 



l]{s,t) 



(1 - 2e)(3 - 2e) 



l X^s,t) l icr{-s)-'l + x^-' 
'4(l + x)2 2 e(l-e) 1 + x 



where x = i/s- The prefactor cr is defined in eq. (|5.6| ), and 

C{s,t) ^ 2(1 - 2e) i-s) (1 + x)lV°'"''''='-''[l]is,t) 

1 x~'^x~'^ ~ (1 ~ ^)~^ 



2icT{-s)-'{l + x) / dx- 
Jo e 

2icr{-s)-' (l + x) C dx 



1-{1 + X)x 



1 



1-{1 + X)x 
Performing the one-dimensional integral we get. 



In 



1 — X 
XX 



(A.6) 



£(ln2(l_x)-ln2(xx))+0(e2) 

(A.7) 



Cis,t) = icri-syl In^x + T^^ 



+ e 



2(- Li3(-x) + C(3) + In X Li2 (-x)) + (In' x + ln(l + x) - 3 In' X " vr' In X 



Using eqs. ( |A.6|) and ( |A.8|) , the desired one-loop integral is 



(A.8) 



r-l-loop r \ 4 



1 J 1 Ixiln^X + T^'^) Xlnx , H 



+ e 



X 



6\e 2 (l + x)2 1 + X ^ 3 
Li3(-X) - C(3) - Inx Li2(-x) + ^ In^ x - ^ In' x ln(l + x) 



+ ^lnp^ 
2 Vl + X 



+ ^((2 + x)ln'x + ^' 



+ 



11 



+ 



+ 0(e'), 



1 X(ln' X + tt') _ xlnx , H 

2 (l + x)2 1 + X 
region (i). 



(A.9) 



The three kinematic regions for the four-point amplitude are 

(i) u > and s,t < 0, for which x > 0; 

(ii) s > and t, u < 0, for which — 1 < x < 0; 

(iii) t > and s, u < 0, for which x < — 1- 
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The form ( [A.9| ) for the one- loop box integral is appropriate for region (i), where it is manifestly real. 
In regions (ii) and (iii) the integral acquires an imaginary part. A form appropriate for region (ii) is 
obtained by substituting (— s)"*^ — > s~'^(l + e z vr — '^^/2) and In x ^ In |x| + ^ tt, 



^l-loop r \ 4 



[Xl]{s,t) = icrs-H-e){l-e] 



ifl 1 Xln Ixl Xln|xl H 



Li3(-X) - C(3) - In Ixl Li2(-x) + I In^ |xl " I In^ Ixl ln(l + x) 



vr^ 1 
+ yln|x| + -(2 + x)ln^|xl 



11 

2(rT^ + Y'"^ 



1 Xln^ Ixl 



X In Ixl 11 
2(l + x)2 1 + x 3 



+ e 



Xln|x| 
(1 + X)2 

X 

(T+x? 



+ 



1 



i + x 



Li2(l + x) + 2 l'^^ Ixl + In Ixl 



11 / Xln|x| 



3V(1+X)2 1 + X 



region (ii). 

(A.IO) 

A form appropriate for region (iii) can be obtained simply by applying the transformation s ^ t 
(x ^ 1/x) to eq. ( A.IO ), since that transformation maps region (ii) into region (iii). 

We now return to the non-vanishing planar double box integral. Its expression in terms of 



r-l-loop r \ 4 



[Ap] is, 



-tcr 



where the finite remainder is 



1 



X 



18(l + x)2 
1 X 
18(l + x)2 



lnx(ln^X + 7r^) + (x ) 



X 



region (i), 



(A.ll) 



In Ixl (ln2 Ixl - 27r2) + (x - i) 

X 



vr^ — Sivr In^ |x| 



regions (ii) and (iii). 

(A.12) 



(A.13) 

-1 in eqs. ( |A1c| ) and ([Al2|) are all 



In rewriting eq. (A. 2), we used 

We remark that the apparent power-law singularities as x 
spurious, cancelling among the various terms. 

The analytic structure of the non-planar integrals is more intricate than that of the planar ones. 
This is because, for Mandelstam variables s, t, u satisfying s + t + u = 0, there is no Euclidean region 
where the integrals are purely real. (It is possible to find a Euclidean region by relaxing the condition 
s + t + u = 0, but this introduces another dimensionless variable and leads to more cumbersome 
expressions |Q.) Also in contrast to the planar case, the non-planar integrals do not have a uniform 
ie prescription in terms of the dimensionless ratio x- The correct ie prescription can be determined 
by following through the ies appearing in the Feynman propagators. 
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We express the divergent non-planar integral in terms of the one- loop box integral, as in eq. ( |A.ll ), 



1 F^^ 



Here the finite remainder is 



P'!r = ^^-2X(1 + X) 



-3(rf^)-C(3)-In(^)(Li.(^)4)-iln3(- 



X 



+ x 



+ 3x(l + X) ln(l + X) In X - ^ (1 + X)' f- - + 3) ln2(l + x) - [-^ 

2 V X / 2 Vl + X 



+ 7r X 



1 1 



+ 7 +(1 + X)ln(l+X)-Xlnx 



+ i7r 2x(l + x) 



2 3 



21 + X 6 

\\^X' 6 2 
region (i), 



Inx 



+ (1 + x) 



(1 + 



+ 3 In^X 



X 



+ 3 ln(l + x)-l 



(A.15) 



Lis 



X 



1 + X 



C(3) - In 



X 



1 + X 



Li2 



X 



1 + X 



6 



In' 



X 



1 + X 



+ 3x(l+x) ln(l+x) Inlxl 



TT 



^(l + x)'(-^ + 3)ln^(l + x) 



+ ^) ^"^^ Ixl + y + (1 + X) ln(l + x) - X In |xl 



-Fzvr x(l +X) 



In^ 



X 



+ vr^ 



(1 + X)(-^ + 2) ln(l + x) + + 2) In Ixl + 1 U , 



1 + X 

region (ii), 

(A.16) 

and the form for region (iii) may be obtained from that for region (i) through the transformation 
t u (x^ ~l~x)- The region (ii) form for is symmetric under x ^ —1 — X; although the 
identities i22|l 



Li3(-x) - Li3(-l/x) 
Li2(-x) + Li2(-l/x) 



^'1 1 1 3 

— m X m X , 

6 6 
^2 



TT 



1 



■In X . 



x > 0, 
X > 0, 



(A.17) 



are needed to demonstrate this. 



Finally, the integral Xl^^fA^ A^] = X^^IA^ A^] is finite, but non-vanishing, as e — > 0. It is given 
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by 



1 



1 



X 



(47r)4(-s)6l(l + x)2 



Li3(-x)-C(3)-lnx(Li2(-x) 



2- 3 
' 4 



X(ln^ X - TT^) 



Lis (3^) - C(3) + ln(l + x)(Li.(3^) + y) + ^(1 + X) In^d + x) 



+ iln3(l+x) 



+ 7r 
+ zvr 



+ 



+ 



+ 



X 



x(i + x) 



6x 3(1 + x) 2(l + x)2 4 



+ 



ln(l + x) Inx 



2X 



2(1 + X) 



1 + X 

■ X^ 



Li2 

X 



X 



2(1 + X)2 



1 + X 
(In' X + vr' 



ln(l + x) Inx + ^ ln2(l + x) " ^(1 + x) ln(l + x) 



1 3 

X(l + X) ^ 2 



Inx 



1 
2^ 



region (i), 



(A.18) 



1 



1 



X 



7-NPrx2 1,2] 

^4 [\+,\\ (4^)4(_,)6|(l + ^)2 



Li3(-X)-C(3)-In|x| Li2(-x) + 



+ ^(^^(YT^ + ^)ln(l + x)ln|x|-- 
X [t- /1 + X 



12x(l + x) S; 2(1 +x) 



4 



Xln' Ixl 



(1+X)2 
1 



Li2 



X 



ln(l +x)ln|x| - 2 I'^lxl 



1 + X 



ln(l + x)-^ln|xl + ^ 



(X 



-1-x) 



region (ii). 

(A.19) 

Again the form for region (iii) may be obtained from that for region (i) through the transformation 
t ^ (x ^ -1 - x)- 
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